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O ; Abstract 

Q^l I shortly review the present status of the theoretical estimates of e'/e . I consider 

p • a few aspects of the theoretical calculations which may be relevant in understand- 

^ ■ ing the present experimental results. In particular, I discuss the role of final state 

interactions and in general of non-factorizable contributions for the explanation of 
the A J = 1/2 selection rule in kaon decays and e'/e . Lacking reliable lattice 
^ ' calculations, the expansion and phenomenological approaches may help in un- 

' derstanding correlations among theoretical effects and the experimental data. The 

same dynamics which underlies the CP conserving selection rule may drive e'/e in 
the range of the recent experimental measurements. 
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Figure 1: The 1-cr results of the NA31 and E731 Collaborations (early 90's) are 
shown by the gray horizontal bands. The old Miinchen, Roma and Trieste theoretical 
predictions for e' je are depicted by the vertical bars with their central values. For 
comparison, the VSA estimate is shown using two renormalization schemes. 

The results announced during the last year by the KTeV 0) and NA48 0) 
collaborations have marked a great experimental achievement, establishing 35 years 
after the discovery of CP violation in the neutral kaon system 0) the existence of a 
much smaller violation acting directly in the decays. 

While the Standard Model (SM) of strong and electroweak interactions 
provides an economical and elegant understanding of indirect [e] and direct (e') CP 
violation in term of a single phase, the detailed calculation of the size of these effects 
implies mastering strong interactions at a scale where perturbative methods break 
down. In addition, CP violation in ^ tttt decays is the result of a destructive 
interference between two sets of contributions, which may inflate up to an order of 
magnitude the uncertainties on the hadronic matrix elements of the effective four- 
quark operators. All that makes predicting e' je a complex and subtle theoretical 
challenge @) . 

The status of the theoretical predictions and experimental data available 
before the KTeV announcement in February 1999 is summarized in Fig. |l]. 

The gray horizontal bands show the one-sigma average of the old (early 
90's) NA31 (CERN) and E731 (Fermilab) results. The vertical lines show the ranges 
of the published theoretical •predictions (before February 1999), identified with the 
cities where most of the group members reside. The range of the naive Vacuum 
Saturation Approximation (VSA) is shown for comparison. 
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Figure 2: The latest theoretical calculations of e' j e are compared with the combined 
l-a average of the NA31, E731, KTeV and NA48 results {e'/e = 19.2 ±4.6 x 10"^), 
depicted by the gray horizontal band (the error is inflated according to the Particle 
Data Group procedure when averaging over data with substantially different central 
values) . 



The experimental and theoretical scenarios have changed substantially af- 
ter the first KTeV data and the subsequent NA48 results. Fig. ^ shows the present 
experimental world average for e'/e compared with the revised or new theoretical 
calculations that appeared during the last year. 

Notwithstanding the complexity of the problem, all theoretical calculations 
show a remarkable overall agreement, most of them pointing to a no n- vanishing 
positive effect in the SM (which is by itself far from trivial). 

On the other hand, if we focus our attention on the central values, many 
of the predictions prefer the 10"'^ regime, whereas only a few of them stand above 
10~^. Is this just a "noise" in the theoretical calculations? 

The answer is no. Without entering the details of the various estimates, 
it is possible to explain most of the abovementioned difference in terms of a single 
effect: the different size of the hadronic matrix element of the gluonic penguin Qq 
obtained in the various approaches. 

While some of the calculations, as the early Miinchen and Rome predic- 
tions, assume for (vrvrlQel-^) values in the neighboroud of the leading result 
(naive factorization), other approaches, among which the Trieste and Dortmund 
calculations, find a substantial enhancement of this matrix element with respect to 
the simple factorization result. The bulk of such an effect is actually a global en- 



hancement of the / = components of the K ^ tttt amphtudes, which affects both 
current-current and penguin operators, and it can be clearly understood in terms of 
chiral dynamics (final-state interactions). 

As a matter of fact, one should in general expect an enhancement of 
e'/e with respect to the naive VSA due to final-state interactions (FSI). As Fermi 
first argued in potential scattering the isospin / = two-body states feel an 
attractive interaction, of a sign opposite to that of the 1 = 2 components thus af- 
fecting the size of the corresponding amplitudes. This feature is at the root of the 
enhancement of the 1 = amplitude over the 1 = 2 one and of the corresponding 
enhancement of e'/e beyond factorization. 

The question is how to make of a qualitative statement a quantitative one. 
Recent dispersive studies of the K ^ nn amplitudes ID make use of the Omnes- 
Mushkelishvili representation 0) 

M{s + le) = P{s) exp (- T — ds'] , (1) 

\ TT J 4ml s' — s — ie j 

in order to resum FSI effects from the knowledge of the vrvr rescattering phase in 
the elastic regime (s < 1 GeV^). -P(s) is a polinomial function of s which is related 
to the factorized amplitude. A solution of the above dispersive relation for the ylo,2 
amplitudes can be written as 

Aj{s)=J^{s-ml)'Ri{s)e'^^^\ (2) 

where A is the derivative of the amplitude at the subtraction point s = m^. The 
coefficient 'K. represents the rescaling effect related to the FSI. By replacing Aj 
with the value given by LO chiral perturbation theory, Pich and Pallante found 
7^(m|)o,2 — 1-4, 0.9 thus confirming via a (partial) resummation of FSI the en- 
hancement of the J = amplitudes, together with a mild depletion of the 1 = 2 
components. 

The numerical significance of these results has been questioned ID on the 
basis that the precise size of the effect depends on the boundary conditions of the 
factorized amplitude which are not unambiguously known, due to the neglect of 
higher order chiral corrections. On the other hand, one should keep in mind that 
the choice of a low subtraction scale minimizes the effect of chiral corrections and 
makes the result of the dispersive analysis trustworthy. 

At any rate a question remains: is the FSI rescaling on the factorized 
amplitudes, albeit a crucial contribution to the calculation of e' j e , all what we were 
looking for ? 



Given the possibility that common systematic uncertainties may affect the 
calculation of e' je and the A/ =1/2 rule (see for instance the present difficulties in 
calculating on the lattice "penguin contractions" for CP violating as well as for CP 
conserving amplitudes ^) a convincing calculation of e' je must involve at the same 
time a reliable explanation of the A/ = 1/2 selection rule, which is still missing. 
FSI effects alone are not enough to account for the large ratio of the 1 = 0,2 
CP conserving amplitudes. Other sources of large non-factorizable corrections are 
needed, which may affect the determination of e' j e as well. 

The A/ = 1/2 selection rule in K tttt decays is known since 45 years m 
and it states the experimental evidence that kaons are 400 times more likely to decay 
in the / = two-pion state than in the 1 = 2 component. This rule is not justified 
by any general symmetry consideration and, although it is common understanding 
that its explanation must be rooted in the dynamics of strong interactions, there is 
no up to date derivation of this effect from first principle QCD. 

The ratio of J = 2 over 1 = amplitudes appears directly in the definition 

of e'/e : 
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As a consequence, a consistent calculation of e'/ e must also address the determina- 
tion of the CP conserving amplitudes. 

The way we approach the calculation of the hadronic K ^ Tin transitions 
in gauge theories is provided by the Operator Product Expansion which allows 
us to write the relevant amplitudes in terms of the hadronic matrix elements of 
effective AS = 1 four quark operators (at a scale /i) and of the corresponding 
Wilson coefficients, which encode the information about the dynamical degrees of 
freedom heavier than the chosen renormalization scale: 

Has=i = ^V^,V:^J2[^,ifi) + Ty,if^)]Q,i^i) . (4) 

The entries Vij of the 3x3 Cabibbo-Kobayashi-Maskawa matrix describe the flavour 
mixing in the SM and r = —VtdV^l/VudV*^. For fj, < rric {q = u,d,s), the relevant 



quark operators are: 

Ql = (SaM/3)v_A (^/3C^a)v-A 

Q7,9 = I {sd)y_j^ J2q (gg)v±A 

aj VitA 



Current- Current 



Gluon "penguins" (5) 



Electroweak "penguins'" 



Current- current operators are induced by tree-level W-excliange whereas the so- 
called penguin (and "box") diagrams are generated via an electroweak loop. Only 
the latter "feel" all three quark families via the virtual quark exchange and are there- 
fore sensitive to the weak CP phase. Current-current operators control instead the 
CP conserving transitions. This fact suggests already that the connection between 
e' /e and the A/ = 1/2 rule is by no means a straightforward one. 

Using the effective AS* = 1 quark Hamiltonian we can write e' /e as 



e' GpUJ 



2|e|ReAo 



U0--U2 

(jj 



(6) 



where 

no = ^EiZ/i Re(g,)o (1 -^^r,+r,') 

and (Qi) = {7nT\Qi\K). The rescattering phases 60,2 can be extracted from elastic 
TT — vr scattering data and are such that cos5o — 0.8 and cos52 — 1. Given that 
the phase of e, 9^, is approximately 7r/4, as well as 60 — 62, (p = f + ^2 — <^o ~ turns 
out to be consistent with zero. 

Two key ingredients appear in eq. |[ 

1. The isospin breaking — rj — t]' mixing, parametrized by flri+ri', which is 
estimated to give a positive correction to the A2 amplitude of about 15-35%. 
However, it has been recently emphasized El) that NLO chiral corrections may 
make negative and as large as —0.7, thus potentially providing a strong 

enhancement mechanism for e'/ e . Further isospin breaking effects have been 
very recently discussed 



2. The combination of CKM elements ImAf = Im(V^*Vtd), which affects directly 
the size of e' / e and the range of the uncertainty. The determination of Im Aj 



goes through B-physics constraints and e. The latter observable depends on 
the theoretical determination of Bk-, the — hadronic parameter, which 
should be consistently determined within every approach. The theoretical 
uncertainty on Bk affects subtantially the final uncertainty on Im Xt. A better 
determination of the unitarity triangle is expected from the B-factories and 
the hadronic colliders 0). From K-physics, Kl ti^vv gives the cleanest 
"theoretical" determination of Im At, albeit representing a great experimental 
challenge 



A satisfactory approach to the calculation of e' /e should comply with the 
following requirements: 

A: A consistent definition of renormalized operators leading to the correct scheme 
and scale matching with the short-distance perturbative analysis. 

B: A self-contained calculation of all relevant hadronic matrix elements (includ- 
ing Bk). 

C: A simultaneous explanation of the A/ = 1/2 selection rule and e' /e . 

None of the available calculations satisfies all previous requirements. I 
summarize very briefly the various attempts to calculate e' /e which have appeared 
so far and have lead to the estimates shown in Figs. |T] and |^. A simple naive 
approach to the problem is the VSA, which is based on two drastic assumptions: 
the factorization of the four quark operators and the vacuum saturation of the 
completeness of the intermediate states. As an example: 

(7r+7r-|Q6|^°) = 2 (7r-|M75rf|0)(7r+|sM|A'°) -2 (7r+7r-|rfrf|0)(0|s75rf|ir°) 



+ 2 



(0|ss|0) - (0|rfrf|0)j (7r+7r-|s75rf|is:°) 



The VSA does not allow for a consistent matching of the scale and scheme depen- 
dence of the Wilson coefficients (the HV and NDR results are shown in Fig. |l]) 
and it carries a potentially large systematic uncertainty iP. It is best used for LO 
estimates. 

Generalized Factorization represents an attempt to address the issue of 
scale and renormalization scheme dependence in the framework of factorization by 
defining effective Wilson coefficients which absorb the perturbative QCD running 
of the quark operators. The new coefficients, scale and scheme independent at a 
given order in the strong coupling expansion, are matched to the factorized matrix 



elements at the scale /ip which is arbitrarily chosen in the perturbative regime. A 
residual scale dependence remains in the penguin matrix elements via the quark 
mass. Fitting the A/ = 1/2 rule and e'/e requires non-factorizable contributions 
both in the current- current and the penguin matrix elements, parametrized by (in- 
dependent) phenomenological parameters . 

In the Miinchen approach (phenomenological 1/A^) some of the matrix 
elements are obtained by fitting the AI = 1/2 rule at /i = rric = 1.3 GeV. On the 
other hand, the relevant gluonic and electroweak penguin (Qe) and {Q8)2 remain 
undetermined and are taken around their leading 1/A^ values (which implies a scheme 
dependent result). In Fig. ^the HV (left) and NDR (right) results are shown 0). 
The dark range represents the result of gaussian treatment of the input parameters 
compared to flat scannning (complete range). 

In the recent years the Dortmund group has revived and improved the ap- 
proach of Bardeen, Buras and Gerard 0) based on the 1/A^ expansion. Chiral loops 
are regularized via a cutoff and the amplitudes are arranged in a p^"/A^ expansion. 
A particular attention has been given to the matching procedure between the scale 
dependence of the chiral loops and that arising from the short-distance analysis 0). 
The renormalization scheme dependence remains and it is included in the final un- 
certainty. The AI = 1/2 rule is reproduced, but the presence of the quadratic cutoff 
induces a matching scale instability (which is very large for Bk)- The NLO correc- 
tions to (Qq) induce a substantial enhancement of the matrix element (right range 
in Fig. 1^) compared to the leading order result (left). The dark range is drawn for 
central values of m^, Qri+ri', ImAj and Aqcd- 

In the Nambu, Jona-Lasinio (NJL) modelhng of QCD H the Dubna 
group has calculated e'/e including chiral loops up to 0{p^) and the effects of 



scalar, vector and axial-vector resonances. Chiral loops are regularized via the heat- 
kernel method, which leaves unsolved the problem of the renormalization scheme 
dependence. A phenomenological fit of the AI = 1/2 rule implies deviations up to a 
factor two on the calculated (Qe)- The reduced (dark) range in Fig. || corresponds to 
taking the central values of the NLO chiral couplings and varying the short-distance 
parameters. 

In the approach of the Trieste group, based on the Chiral Quark Model 
(xQM) S', all hadronic matrix elements are computed up to 0{p^) in the chiral 
expansion in terms of the three model parameters: the constituent quark mass, the 
quark condensate and the gluon condensate. These parameters are phenomeno- 
logically fixed by fitting the AI = 1/2 rule This step is crucial in order to 
make the model predictive, since there is no a-priori argument for the consistency 



of the matching procedure. As a matter of fact, all computed observables turn 
out to be very weakly dependent on the scale (and the renormalization scheme) in 
a few hundred MeV range around the matching scale, which is taken at 0.8 GeV 
as a compromise between the ranges of validity of perturbation theory and chiral 
expansion. The 1 = matrix elements are strongly enhanced by non-factorizable 
contributions and drive e'/e in the 10~^ regime. The dark (light) ranges in Fig. |^ 
correspond to Gaussian (flat) scan of the input parameters. The bars on the left 
represent the results of ref. EH) which updates the 1997 calculation EID. That on 
the right are a new estimate 0) , similarly based on the xQM hadronic matrix ele- 
ments, in which however e'/e is estimated in a novel way by including the explicit 
computation of e in the ratio as opposed to the usual procedure of taking its value 
from the experiments. This approach has the advantage of being independent from 
the determination of the CKM parameters Im At and of showing more directly the 
dependence on the long-distance parameter as determined within the model. 

Lattice regularization of QCD is the consistent approach to the problem. 
On the other hand, there are presently important numerical and theoretical lim- 
itations, like the quenching approximation and chiral symmetry, which may sub- 
stantially affect the calculation of the weak matrix elements. In addition, chiral 
perturbation theory is needed in order to obtain K ^ amplitudes from the 
computed K ^ tt transitions. As summarized by Ciuchini at this conference 1^) 
lattice cannot provide us at present with reliable calculations of the 1 = penguin 
operators relevant to e'/e , as well as of the / = components of the hadronic 
matrix elements of the current- current operators (penguin contractions), which are 
relevant to the A/ = 1/2 rule. This is due to large renormalization uncertainties, 
partly related to the breaking of chiral symmetry on the lattice. In this respect, very 
promising is the Domain Wall Fermion approach 0) which allows us to decouple the 
chiral symmetry from the continuum limit In the recent Roma re-evaluation 

of e'/e (Qq) is taken at the VSA value with a 100% uncertainty HP. The result is 
therefore scheme dependent (the HV and NDR results are shown in Fig. |^). The 
dark (light) ranges correspond to Gaussian (fiat) scan of the input parameters. 

The A/ = 1/2 rule and Bk have been studied in the NJL framework and 
1/A^ expansion by Bijnens and Prades 0) showing an impressive scale stability 
when including vector and axial-vector resonances. The same authors have recently 
produced a calculation of e'/e at the NLO in 1/A^ ilD. The calculation is done 
in the chiral limit and it is eventually corrected by estimating the largest SU{3) 
breaking effects. Particular attention is devoted to the matching between long- and 
short-distance components by use of the X-boson method ^3) . The couplings of 
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Figure 3: Predicting e' je: a (Penguin) Comparative Anatomy of the Miinchen (dark 
gray) and Trieste (light gray) resuhs in Fig. |1| (in units of 10~^). 



the X-bosons are computed within the ENJL model which improves the high-energy 
behavior. The A/ = 1/2 rule is reproduced and the computed amphtudes show a 
satisfactory renormalization scale and scheme stability. A sizeable enhancement of 
the Qe matrix element is found which brings the central value of e' j e at the level of 
3 X 10-3. 

Other attempts to reproduce the measured e' /e using the linear cr-model, 
which include the effect of a scalar resonance with m„ ~ 900 MeV, obtain the needed 
enhancement of (Qg) ■ However, it is not possible to reproduce simultaneously 
the experimental values of e' /e and of the CP conserving K irir amplitudes. 

Studies on the matching between long- and short- distances in large N 
QCD, with the calculation of the Qr penguin matrix element and of Bk at the 
NLO in the 1/A^ expansion have been presented in ref. However, a complete 
calculation of the K ^ tttt matrix elements relevant to e'/e is not available yet. 

Without entering into the details of the various calculations I wish to il- 
lustrate with a simple exercise the crucial role of final-state interactions (in general 
of non-factorizable contributions) for e'/e and the A J =1/2 selection rule. In order 
to do that I focus on two semi-phenomenological approaches. 

A commonly used way of comparing the estimates of hadronic matrix el- 
ements in different approaches is via the so-called B factors which represent the 
ratio of the model matrix elements to the corresponding VSA values. However, care 
must be taken in the comparison of different models due to the scale dependence of 
the B's and the values used by different groups for the parameters that enter the 



3.5 



> 
O 



00 

< 
o 



2.5 



(N 
II 



< 1.5 



(1) 



(2) 



0.5 



(3,4) 



(5) 



(6) 

• (8) 



(7) J 



1.5 2 2.5 

A(I=0) X 10^7 (GeV) 



3.5 



Figure 4: Anatomy of the AJ = 1/2 rule in the xQM. See the text for explanations. 
The cross-hairs indicate the experimental point. 



VSA expressions. An alternative pictorial and synthetic way of analyzing different 
outcomes for e'/e is shown in Fig. where a "comparative anatomy" of the old 
Trieste and Miinchen predictions is presented. 

Prom the inspection of the various contributions it is apparent that the 
different outcome on the central value of e' / e is almost entirely due to the difference 
in the Qq component. 

In the Miinich approach 113) the A J = 1/2 rule is used in order to determine 
phenomenologically the matrix elements of Qi,2 and, via operatorial relations, some 
of the matrix elements of the left-handed penguins. The approach does not allow for 
a phenomenological determination of the matrix elements of the penguin operators 
which are most relevant for e'/ e , namely the gluonic penguin Qq and the electroweak 
penguin Qg. These matrix elements are taken around their leading 1/A^ values 
(factorization). 

In the semi-phenomenological approach of the Trieste group the size of the 
effects on the J = 0, 2 amplitudes is controlled by the phenomenological embedding 
of the A/ = 1/2 selection rule which determines the ranges of the model paremeters: 
the constituent quark mass, the quark and the gluon condensates. In terms of these 
parameters all matrix elements are computed. 

Fig. ^ shows an anatomy of the xQM contributions which lead to the 
experimental value of the A/ = 1/2 selection rule for central values of the input 
parameters. Point (1) represents the result obtained by neglecting QCD and taking 
the factorized matrix element for the tree-level operator Q2, which is the leading 



electroweak contribution. The ratio A0/A2 is thus found equal to ^/2: by far off the 
experimental point (8). Step (2) includes the effects of perturbative QCD renormal- 
ization on the operators (5i,2 ■ Step (3) shows the effect of including the gluonic 
penguin operators . Electroweak penguins are numerically negligeable in the 
CP conserving amplitudes and are responsible for the very small shift in the A2 
direction. Therefore, perturbative QCD and factorization lead us from (1) to (4). 

Non-factorizable gluon-condensate corrections, a crucial model dependent 
effect entering at the leading order in the chiral expansion, produce a substantial 
reduction of the A2 amplitude (5), as it was first observed by Pich and de Rafael . 
Moving the analysis to O(p^), the chiral loop corrections, computed on the LO chiral 
lagrangian via dimensional regularization and minimal subtraction, lead us from 
(5) to (6), while the finite parts of the NLO counterterms calculated in the xQM 
approach lead to the point (7). Finally, step (8) represents the inclusion of vr- 77-77' 
isospin breaking effects . 

This model dependent anatomy shows the relevance of non-factorizable 
contributions and higher-order chiral corrections. The suggestion that chiral dy- 
namics may be relevant to the understanding of the A/ = 1/2 selection rule goes 
back to the work of Bardeen, Buras and Gerard — ^ in the framework with a 



cutoff regularization. A pattern similar to that shown in Fig. ^ for the chiral loop 
corrections to Aq and A2 was previously obtained in a NLO chiral lagrangian anal- 
ysis, using dimensional regularization, by Missimer, Kambor and Wyler The 
Trieste group has extended their calculation to include the NLO contributions to 
the electroweak penguin matrix elements S HI) . 

Fig. 1^ shows the contributions to e'/ e of the various penguin operators, pro- 
viding us with a finer anatomy of the NLO chiral corrections. It is clear that chiral- 
loop dynamics plays a subleading role in the electroweak penguin sector (Qs-io) 
while enhancing by 60% the gluonic penguin (/ = 0) matrix elements. The NLO 
enhancement of the Qq matrix element is what drives e'/e in the xQM to the 10~^ 
ballpark. 

As a consequence, the xQM analysis shows that the same dynamics that is 
relevant to the reproduction of the CP conserving Aq amplitude (Fig. ^) is at work 
in the CP violating sector, albeit with a reduced strenght. 

In order to ascertain whether these model features represent real QCD 
effects we must wait for future improvements in lattice calculations S 0) . On the 
other hand, indications for such a dynamics follow also from the analitic properties 
of the K ^ nn amplitudes, as the dispersive analyses show ^) . It is important to 
notice however that the size of the effect so derived is generally not enough to fully 
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Figure 5: Anatomy of e'/e (in units of 10~^) in the xQM approach. In black the 
LO results (which includes the non-factorizable gluonic corrections), in half-tone the 
effect of the inclusion of chiral-loop corrections and in light gray the complete 0{p'^) 
estimate. 



account for the A J = 1/2 rule. Other non-factorizable contributions are needed 
to further enhance the CP conserving / = amplitude, and to reduce the large 
1 = 2 amplitude obtained from perturbative QCD and factorization. In the xQM 
approach, for instance, the fit of the AI = 1/2 rule is due to the interplay of FSI 
(at the NLO) and non-factorizable soft gluonic corrections (at LO in the chiral 
expansion) . 

The idea of a connection between the AI = 1/2 selection rule and e'/e goes 



back a long way although at the GeV scale, where we can trust perturbative 
QCD, penguins are far from providing the dominant contribution to the CP con- 
serving amplitudes. 

In summary, those semi-phenomenological approaches which reproduce the 
AI = 1/2 selection rule in K —>■ decays, generally agree in the pattern and size 
of the 1 = 2 hadronic matrix elements with the existing lattice calculations. On the 
other hand, the AI =1/2 rule forces upon us large deviations from the naive VSA 
for the / = amplitudes: 5— factors of 0(10) are required for {Qi^2)q- Here is were 
lattice calculations suffer from large sistematic uncertainties. 

In the Trieste and Dortmund calculations, which reproduce the CP con- 
serving K ^ Tin amplitudes, non-factorizable effects (mainly final-state interactions) 



(2) _ 9 
8 ~ ^• 

Similar indications stem from recent l/N^^^ and dispersive approaches. The 



enhance the hadronic matrix elements of the gluonic penguins, and give Bq/B, 



direct calculation of K ^ vrvr amplitudes and unquenching are needed in the lattice 
calculations in order to account for final state interactions. Further progress in this 
direction is awaited. 
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